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We consider pure isocurvature cold dark matter models in the case of open and closed universes. 
We allow for a large spectral tilt and scan the 6-dimensional parameter space for the best fit to the 
COBE, Boomerang, and Maxima-1 data. Taking into account constraints from large-scale structure 
and big bang nucleosynthesis, we find a best fit with = 121, which is to be compared to = 44 
of a flat adiabatic reference model. Hence the current data strongly disfavor pure isocurvature 
perturbations. 



I. INTRODUCTION 

The recent measurements of the cosmic microwave 
background (CMB) temperature fluctuations by the 
Boomerang and Maxima-1 0,^ balloon experiments 
and the DASI interferometer Q have widely been re- 
garded as indicating that we live in a = 1 universe. 
This is so because the first acoustic peak is found at the 
multipole (. ~ 200, implying a flat universe. The firmness 
of such a conclusion is, however, based on certain tacit as- 
sumptions. In particular, when fitting the acoustic peak 
positions, one often assumes that the primordial pertur- 
bations are adiabatic and that the spectrum is nearly 
scale invariant. 

If perturbations are adiabatic, the relative abundances 
of particle species are equal to their thermal equilibrium 
values. This is the case in the simplest, one-field infla- 
tion models but it is not a generic feature of inflation. 
More generally, perturbations can be either adiabatic or 
nonadiabatic; the latter would be perturbations in the 
particle number densities, or entropy perturbations, and 
are called isocurvature perturbations. 

Because no generally accepted theory of inflation ex- 
ists, it is natural to consider both adiabatic and isocur- 
vature perturbations as being equally probable. This is 
the generic situation when more than one field is excited 
during inflation, such as is the case in double inflation 
Q or in the minimally supersymmetric standard model 
with flat directions 0. One should also note that in 
the pre-big-bang scenario, which has been proposed as 
an alternative to the inflationary universe, pre-big-bang 
axion field fluctuations give rise to an isocurvature per- 
turbation spectrum |^]. Purely isocurvature f2 = 1 per- 
turbations are, however, not consistent [p|-pl| with the 



observational data, but an admixture of (uncorrelated 
or correlated) adiabatic and isocurvature perturbations 
cannot be ruled out Jll|-p^. However, if we do not insist 
on a flat universe, the situation could be different. 

Recently, it was pointed out that in the general 
(Gaussian) case the scalar power spectrum is a 5 x 5 
matrix Pij(k) = (Ai(k)Aj(— k)), where i, j label one adi- 
abatic and four isocurvature modes [cold dark matter 
(CDM), baryon, neutrino density, and neutrino velocity] 
and their correlations. Here we shall focus on a purely 
isocurvature primordial perturbation in the CDM which 
has the power spectrum 



Ps(k) = B (^\/k^ - K 



(1) 



where niso is the spectral index and K — ~Hq{1 — fl) 
is the curvature. Since in curved space the Laplacian 
has eigenvalues — K instead of the /c^ of the flat case, 
the spectrum (0) is the simplest generalization of the flat 
space spectrum 

In the flat, = 1 case, definition (|l|) gives the 
power law Ps{k) oc k""°~^, which is a natural form for 
the power spectrum, and approximates well the spec- 
trum produced by typical inflation models with isocur- 
vature perturbations in the region of interest. The scale- 
invariant spectrum has njso = 1. In principle, riiso could 
well depend on fc; here we shall assume that it is a con- 
stant (or varies very little) over the range of interest. In 
open and closed models the spatial curvature introduces 
a length scale and one expects this to be reflected in the 
form of the power spectrum. It is not obvious what would 
be the most natural modification to the power law for 
isocurvature models in curved space. This question has 
been studied only for specific models in the adiabatic case 
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1^,0. Thus we stress that we are using a phenomeno- 
logical power-law spectrum, which does not necessarily 
follow from any particular inflation model. We shall re- 
turn to this point later in this paper. 

After the clear detection of the acoustic peak around 
I ~ 200 it became evident that the adiabatic models fit 
weh to the data ||l|,|,|j|Jl|,|l|l . However, this should not 
be taken as a proof that all pure isocurvature models are 
ruled out. Some unconventional combination of cosmo- 
logical parameters, e.g., 7^ 1 and a spectrum with a 
large tilt, could at least in principle give an equally good 
fit as do the adiabatic models. 

Pure isocurvature models have two well-recognized 
problems: excess power at low multipoles and a peak 
structure that is roughly speaking out of phase by 7r/2 
when compared to the adiabatic one pof . Since the angu- 
lar power in the low multipole region was measured quite 
firmly by the Cosmic Background Explorer (COBE), ^ 
fitting forces the overall normalization constant in pure 
isocurvature models to be smaller than in the adiabatic 
case, which leads to too little power at higher multipoles. 
The easiest and perhaps the only way to compensate for 
this is to introduce a large spectral tilt. Moreover, since 
flat adiabatic models fit the observed peak at £ ~ 200 
well, it is obvious that the I ~ 200 peak falls between the 
first and second peaks of any flat isocurvature model. 
Accordingly, in our earlier study [0, the best-fit flat 
isocurvature model was found to have a large — 116 
for 30 data points and 6 parameters whereas the best 
adiabatic model had = 22. 

Thus we have two possibilities for a better isocurva- 
ture model. The first is to lower the total energy density 
parameter so much that the position of the first isocur- 
vature peak fits to the observed peak at £ ~ 200, which 
means that we have to allow for an open universe (11 < 1). 
The other possibility is to increase the total energy den- 
sity parameter so much that the position of the second 
isocurvature peak fits the i ~ 200 peak j2^, implying a 
closed universe (17 > 1). In this case the first isocurva- 
ture peak at £ ~ 60 . . . 100 should effectively disappear. 
In fact, a large spectral tilt would have precisely this ef- 
fect since it would decrease the relative power at low t. 

The purpose of the present paper is to study these pos- 
sibilities systematically to find out if CDM isocurvature 
models are indeed completely ruled out by the presently 
available CMB data. 



II. METHODS AND RESULTS 

In order to compare the isocurvature models 
with adiabatic ones we choose one representative 
well-fitted adiabatic model (riadi, l^m, l^A, i^b, i^c, ''') = 
(0.98,0.38,0.62,0.021,0.13,0); cf. 0. Using the same 
data sets and algorithm as for isocurvature models, 
we get = 44 for this adiabatic "reference" model. 
Fig. ||(b) confirms that this model fits well both the low 



I part of the angular power spectrum and the acoustic 
peaks. 

Our starting point for analyzing isocurvature models 
is a large grid with the following free parameters: 

• riiso = 1.00 . . . 7.00 (60 values) 

• f7„ 0.06 . . . 2.31 (16 values) 

• = -1.00 . . . 1.10 (14 values) 

• wfc = 0.001 . . .0.100 (10 values) 

• ijJa = 0.01 ... 1.60 (15 values), 

where 17m is the total matter density, 17a is the vacuum 
energy density, tjf, = /i^l7f, is the baryon density, and Wc = 
/i^l7c is the CDM density. The sixth free parameter is the 
overall normalization factor B of Eq. (P. The Hubble 
constant h is not a free parameter, since h^flm = = 
bjfj -I- ujf.- We use a top-hat prior h = 0.45 . . .0.90 and 
assume r = for the optical depth due to reionization. 
The angular power spectrum of all the models in the 
grid was calculated by CAMB assuming isocurvature 
CDM initial conditions. 

We use the x^ method to compare models and data, 
because it allows us to quickly search a large parameter 
space. This method is approximate and we do not 
attempt precise estimates for cosmological parameters or 
confidence levels. As will be seen, the conclusion is clear 
enough in ruling out the isocurvature models so that it is 
not necessary to go to a full maximum likelihood analysis 

Using the latest Boomerang data [|1|, together with 
Maxima- 1 |^ and COBE data Q we calculate x^ for 
each model. The resulting best-x^ contours in the 
(17m, 17a) plane are presented in Fig. |l| by gray levels. 
The best-fit model turns out to have x^ — 80 with 
(niso,17m,17A,W6,Wc) = (2.00,2.11,-1.00,0.020,1.40). 
From Fig. [^(a) we see that the best-fit isocurvature mod- 
els lie along two bands in the (17m, 17a) plane, the left 
band corresponding to open universes, and the right cor- 
responding to closed universes. In the best- fit models the 
spectral index falls in the range niso = 2 ... 3. 

A detailed examination of the various pure isocurva- 
ture models allows us to conclude that the main prob- 
lems are the spacings of the higher acoustic peaks and 
the slope in the (low £) Sachs- Wolfe region. COBE mea- 
sured a close-to-flat Q spectrum, but the isocurvature 
models have a signiflcant positive slope arising from the 
large primordial blue spectral tilt needed to get enough 
power at higher multipoles. 

In the best-fit open models the prominent peak in the 
CMB data is fitted by the first acoustic peak of the isocur- 
vature model. Fig. |l|(a) shows that in the best-fit open 
region the first peak lies in the range 150 < ^ < 230. 
Since the data do not show a high second peak, these 
models need a small baryon density Ub to boost up the 
first peak and suppress the second peak. (In the adia- 
batic case, adding more baryons enhances odd acoustic 
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FIG. 1. The hest-x^ contours on the {0,m,0,A) plane. The 
best fit, which has = 80, is indicated by an asterisk (*) 
near to the lower right corner. The contours for deviation 
fi-om the best fit are as follows: white Ax^ < 10; light gray 
10 < Ax^ < 40; medium gray 40 < Ax^ < 100; and dark 
gray Ax^ > 100. (a) Dashed lines show the position {£) of 
the first acoustic peak and solid lines the second peak, (b) 
Solid lines give the values of agO,'^^, and the dotted area is 
that allowed by the LSS constraint 0.43 < asQ.°;,^^ < 0.70. 
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FIG. 2. (a) As Fig. g(a) but now with the LSS constraint 
0.43 < cTgnJl;*'^ < 0.70. The best fit marked by an asterisk has 
X^ ~ 103. The contours for deviation from the best fit are as 
follows: white Ax^ < 35; light gray 35 < Ax'^ < 140; medium 
gray 140 < Ax^ < 350; and dark gray Ax^ > 350. The upper 
left corner corresponds to the closed models where the second 
acoustic peak fits the prominent peak in the Ce data, (b) The 
best-fit physical region using the fine grid. The solid contours 
show the baryon density cji,. The best-fit model has x^ = 121 
and the gray levels are as follows: white Ax^ < 6; light gray 
6 < Ax^ < 30, medium gray 30 < Ax^ < 60, and dark gray 
Ax^ > 60. 



peaks over even pX^ , but in the isocurvature case increas- 
ing u)b boosts even peaks.) Actually, all the best-fit open 
models have a baryon density oiojb = 0.001, which is the 
smallest value in the grid. However, even assuming such 
an unphysically low baryon density as 0.0005 only gives 
about half of the power needed to fit the first peak, so 
not scanning below ojb < 0.001 seems justified. 

In the best-fit closed models the i ~ 200 peak in 
the CMB data is fitted by the second isocurvature 
peak, which lies, according to Fig. |l|(a), in the range 



225 < £ < 265. As one might expect (see, e.g., [|5) for 
an adiabatic analogy) , now the ratio of the £ ~ 200 peak 
to the higher multipole C'iS in the data fixes uji, near the 
value 0.02 in the whole best-fit band. In contrast one ob- 
tains almost no restriction for lUc- This is consistent with 
Fig. |l], where fim can be seen to be able to take almost 
any value, which is then compensated by Ha to produce 
the correct peak position. 

According to Fig. |^(a) the best isocurvature model 
{x^ = 80) does badly with the COBE region as well as 
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FIG. 3. Angular power spectra for different models along 
with COBE (o), Boomerang (•) , and Maxima-1 (o) data, (a) 
Best-fit isocurvature model of Fig. ^ (solid line) and best-fit 
open model with LSS constraint (dashed line), (b) Best phys- 
ical isocurvature fit from the fine grid (solid line) and the adi- 
abatic reference model (dashed line) . Note that up to £ = 25 
the I axis is logarithmic. 



after the prominent peak. This peak is fitted quite well 
by the second acoustic peak while the first acoustic peak 
appears as a small shoulder around I ~ 80. 

The considerations so far rely on the CMB data only. 
However, as is well known, when discussing isocurva- 
ture models it is essential to include also the large- 
scale structure (LSS) data. As we will see, rough mea- 
sures are already very effective in constraining the mod- 
els. Therefore we make use of the the amplitude of the 
rms mass fluctuations in an 8/i~^ Mpc sphere only, de- 
noted as (Ts, which the LSS data restricts to the range 
0.43 < crsf^m^^ < 0.70 1^. The contours of a^^^^^ are 
shown in Fig. |l|(b). Apart from the upper left corner of 
the (f7„j,riA) plane, the best-fit closed models appear 



to give a far too large a%^^^ > 1.5. This is natural, 
since we need a large njso to do away with the first peak 
("the isocurvature shoulder") at £ ~ 60 . . . 100 and to get 
enough power at higher multipoles. A large njso evidently 
leads to a large as- To compensate for this, one would 
require a small ^Im- We have checked that the smaller 
flm we have, the larger niso is allowed for by the LSS con- 
straint. In particular, the upper left corner closed models 
in Fig. |l|b obey the LSS constraint, although they have 
a rather large spectral index njso ~ 3.1. 

On the other hand, the best-fit open models tend to 
have a slightly too small as^m^^- These models have a 
relatively small riiso S 2.1, for the following reasons. (1) 
Since these models fit the first isocurvature peak to the 
£ ~ 200 peak in the data, they do not need a large njso 
to eliminate this first peak. (2) The smaller scales do not 
need as large a boost from niso, since power is provided 
by the second peak where the data requires it. Because 
of this smaller riiso these models fit the COBE region 
better. 

We have repeated the analysis of minimizing but 
now with the LSS constraint. As one might expect, this 
eliminates most of the best-fit closed models, leaving only 
those with a small il^ and a large J^a; see the upper left 
corner of Fig. ^(a). The reason for this shifting of the 
best-fit closed-model region to the opposite corner in the 
(^rm^A) plane is easy to understand. Large riiso leads 
to a large as, and hence the prior 0.43 < asi^^^^ < 0.70 
requires flm to be small, which in turn implies a large 
f2A in order to adjust the peak position. 

After imposing the LSS constraint, the best-fit model 
is no longer a closed one but an open model at the cor- 
ner of the parameter space with lji, = 0.001 and Oa = 
— 1.00. This fit has — 103 and (niso, f^m, ^A, ijJb, ijJc) = 
(2.05,0.71,-1.00,0.001,0.16). Fig. |(a) shows that the 
first acoustic peak at £ ~ 170 is too low to fit the data. It 
is clear that the fit would further improve if one allowed 
for even smaller LUh and However, such a small Ub 
is in clear conflict with big bang nucleosynthesis (BBN) . 
There is some debate in the BBN community on 
how small an Uh could be acceptable. After imposing 
a very conservative lower limit, ivb > 0.003, our best-fit 
open model is already significantly worse than the best-fit 
closed models. Moreover, the best-fit open models have 
a very small, even a negative, This region of the 

(firm ^ a) plane is disfavored by the observed supernova 
redshift-distance relationship p^ . 

Thus we conclude that the best candidates for pure 
isocurvature models are the remaining best-fit closed 
models. These models satisfy the LSS constraint and 
have an acceptable ujh. They lie in the region of 
small flm and large Oa. We scanned this region with 
a finer grid. The resulting best-v^ contours in the 
(rim,riA) plane are shown in Fig. ||(b) along with the 
baryon density of these models. The best "physically 
acceptable" isocurvature fit has (riiso, f^m, f^A, wt, o^c) — 
(2.80,0.12,0.97,0.015,0.074). The fit remains very bad, 
however, with = 121 for 40 data points and 6 pa- 
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rameters, to be compared to = 44 of the flat adia- 
batic reference model. Because of the high of the best 
fit, it is unnecessary to consider the LSS spectrum in a 
more detailed way. The badness of the fit is mainly due 
to the COBE and Boomerang data; see Fig. ^(b). The 
COBE contribution to is 2.4 per COBE data point, 
the Boomerang contribution is 4.2 per data point, while 
the Maxima contribution remains at 1.7. The slope of the 
best-fit model is the reason for the poor fit to COBE, and 
although the prominent peak in the data is fitted quite 
well, the "flat adiabatic" peak structure of the second 
and third peaks in the Boomerang data leads to a con- 
flict with the isocurvature peak structure. 

As mentioned earlier, the power-law form for the power 
spectrum is not necessarily the most natural one in open 
and closed models due to the effect of spatial curvature. 
The curvature scale in the models studied is compara- 
ble to the Hubble scale, or larger. Thus its effect is ex- 
pected to be reflected in the COBE region of the power 
spectrum, but not in the Boomerang/Maxima region. 
To assess the significance of this problem, we repeated 
our analysis without the COBE data points. The re- 
sults remained essentially unchanged. Without the 8 
COBE points we got = 70 for the best-fit model, 
X^ = 91 for the best-fit with LSS constraint, — 89 for 
the best physically acceptable fit from the refined grid, 
and x^ — 40 for the adiabatic reference model. Hence 
the Boomerang data alone are sufficient to rule out pure 
isocurvature models and our conclusions do not depend 
on the question of the effect of spatial curvature on the 
power spectrum. 

Actually, since the main discriminant is the relative po- 
sitions of the three peaks in the Boomerang data, which 
show an "adiabatic" instead of an "isocurvature" pat- 
tern, our conclusion should be independent of the shape 
of the primordial power spectrum as long as the observed 
peaks are indeed due to acoustic oscillations and do not 
represent features of the primordial power spectrum it- 
self. 



III. SUMMARY 

We have surveyed a large space of parameters for 
pure isocurvature models, and allowed for both open and 
closed universes, to find out whether there are any pure 
isocurvature models that fit the current CMB data better 
than or at least equally as well as the flat adiabatic model. 
There are none. We conclude that, even if one ignores 
the high-2 supernova data, pure isocurvature CDM mod- 
els, including the ones with a heavily tilted spectrum, are 
completely ruled out by the present CMB and LSS data. 
Incidentally, the isocurvature models do not do too badly 
with the Maxima- 1 data. The main CMB problems are 
with the COBE and the Boomerang data. To have suffi- 
cient smaller-scale power, and to suppress the first peak 
and boost the second peak in the closed models, a large 



blue tilt is needed. This leads to a slope in the Sachs- 
Wolfe region and reduces the largest-scale power below 
the level observed by COBE. The most significant prob- 
lem, however, is with the Boomerang data. Boomerang 
shows a second and a third peak with a spacing that cor- 
responds to a flat universe, whereas the position of the 
flrst peak in the data cannot be fltted by flat isocurvature 
models. 
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